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It is shown how the gravity-gradient effect may be utilized to design a 
long-lived, earth-pointing satellite attitude control system which requires 
no fuel supplies, attitude sensors or active control equipment. This two- 
body system is provided with a magnetic hysteresis damper which effectively 
damps out oscillations (librations) about the local vertical. The long rods, 
which must be extended in space from coiled up metal tapes, provide the 
required large moments of inertia and possess adequate rigidity and suffi- 
cient strength to endure the rigors of the extension process. The system is 
compatible with the requirements of multiple satellite launchings from a 
single last-stage vehicle. Analysis indicates that the gravitational torques 
are sufficient to keep the disturbing effects of solar radiation pressure, re- 
sidual magnetic dipole moments, orbit eccentricity, rod curvature, eddy 
currents, and meteorite impacts within tolerable limits. It is believed that 
the high-performance, earth-pointing system described and analyzed in 
this paper represents an essential step in the development of high-capacity 
communications satellites requiring long life. 
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I. INTRODUCTION 

An earth-pointing attitude control system offers many advantages 
for a commercial satellite repeater. By directing the satellite's radiated 
power to just cover the earth, the satellite's size and weight can be 
minimized. For example, at a 6000-nm altitude the theoretical gain of 
an earth-covering conical beam is 14.5 db; however, allowance must be 
made for inaccuracies of the earth-pointing system and the gain that 
can be achieved from a practical antenna. Conservative estimates have 
shown that the achievable antenna gain is at least 10 db higher than 
with a Telstar-type isotropic antenna.* Hence, with an earth-pointing 
antenna, the power required from the satellite transmitter is only one- 
tenth that required with an isotropic antenna. This reduction in power 
makes the size and weight of communications satellites of high capacity 
(e.g., two TV channels or 600 two-way voice channels continuously 
operating) compatible with existing launch vehicles for orbits of interest. 

In this paper we will describe a passive gravitational attitude control 
system (hereafter called PGAC) which provides a particularly attractive 
way to maintain a satellite axis pointing towards the earth. This system 
should have an extremely long life since it is entirely passive and re- 
quires no power and no active controls or attitude sensors. The system 
has been designed to be compatible with the launching of several satel- 
lites from a single launch vehicle. The importance of this feature becomes 

* While the first Telstar satellite satisfied the objectives for a communications 
experiment, the performance was about 6 db below Bell System objectives. 1 For 
a higher-altitude commercial satellite, the additional 10 db would be considered 
essential to assist in meeting systems margins. 
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apparent when one considers the cost and time required to place perhaps 
20 or more satellites into orbit with existing launch vehicles for a me- 
dium-altitude satellite system. Unpublished studies at Bell Telephone 
Laboratories have indicated that three satellites of the capacity men- 
tioned and employing PGAC can be launched by a single Atlas-Agena 
vehicle in orbits of communications interest. 

The most critical part of any passive earth-pointing system is the 
technique of damping employed to stop tumbling and limit librational 
motions. A unique feature of the PGAC system described herein is the 
employment of magnetic hysteresis damping in conjunction with a 
two-body system that provides large relative motion for damping 
purposes. Magnetic hysteresis damping is quite effective even at the 
slow librational rates (approximately a six-hour period at 6000-nm 
altitude). 

PGAC employs long extensible rods to obtain appropriate moments 
of inertia about the three principal axes. In this respect it is similar to 
other passive systems 2,3 which also require long rods to obtain a suffi- 
ciently large moment of inertia in order for the gravity torque to be 
effective. Another feature of this PGAC system is that a single trigger 
or signal separates each satellite and simultaneously causes the rods to 
extend. This simplicity should enhance reliability of satellite separation 
and rod extension. 

In any passive gravity-gradient orientation system, the satellite is 
stable with either end pointing towards the earth. In the PGAC system 
described here, dual antennas are proposed for each end of the satellite, 
and the appropriate antennas are to be activated by a simple microwave 
switch. Fig. 1 shows the two possible stable positions of the satellite. 
Antenna tests have shown that the extended rods do not substantially 
affect the antenna pattern; the maximum loss due to the rods is about 
1 db. However, it may be possible to avoid this loss by properly orienting 
each satellite initially. This would require precise control of the launching 
vehicle orientation, satellite tumbling rate during ejection, and speed of 
extension of the rods. 

In Section II of this paper, the dynamic principles of PGAC are 
described, and a general description of the system is given. Actually, 
two alternative configurations are described, each of which has its own 
advantages. Vibration analysis of the system is then given in Section 
III to demonstrate the validity of certain rigid body assumptions made 
in the dynamics analysis of the accompanying paper. 4 The stress and 
deformation of the rods due to dynamic loading during the extension 
phase and due to thermal effects are analyzed in Sections IV and V. 
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Fig. 1 — Possible satellite orientations. 



Various spring designs for satellite separation associated with multiple 
launch are described in Section VI. The status of the hardware develop- 
ment and tests on the hysteresis damper unit are reviewed in Section 
VII. Finally, the various disturbing torques which the satellite will 
encounter in space are reviewed in Section VIII. It is shown that the 
PGAC system should remain earth-pointing within a few degrees. 

Typical computer results have disclosed that for a reasonable initial 
tumbling rate of the satellite (1 rpm before rod extension, due to ejection 
from the rocket), the satellite will be earth-pointing within a few degrees 
of the local vertical in about 10 to 15 orbital periods. The description 
and discussion of PGAC in this paper is primarily for a satellite in a 
circular 6000-nm orbit with any inclination. However, with modifica- 
tions of rod lengths and damping and spring constants, PGAC could 
be adapted to either higher or lower orbits. 

The companion paper 4 in this issue covers the basic dynamics analysis 
of PGAC. The analysis includes large angle motion (as would be experi- 
enced by a satellite due to tumbling after ejection from the launch 
vehicle), as well as small librational motion. A complete three-dimen- 
sional analysis of the satellite motion has been formulated, and stability 
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criteria for the system have been determined. Other dynamics analyses 
of passive attitude control systems have been reported in the literature. 5 ' 6 
These analyses either have not included large angle motions or have been 
restricted to the pitch motion only. 

II. DYNAMIC PRINCIPLES AND GENERAL DESCRIPTION OF PGAC 

2.1 Principles 

The fact that an elongated body in orbit around the earth tends to 
line up with the local vertical is well known. 7 Just why this should be 
so is most easily explained by considering a rigid dumbbell with equal 
tip masses. Fig. 2(a) shows a dumbbell, in orbit around the earth, whose 
axis makes an angle 0(0 < 90°) with the local vertical. Since the gravi- 
tational attraction varies inversely as the square of the distance from 
the geocenter, the lower mass A will experience a gravity force F A which 
is slightly larger than the force F a experienced by the upper mass B. 



ORBIT 




BQtip MASS 




SATELLITE 



(C) 



4 



Fig. 2 — (a) Gravity forces acting on a dumbbell in orbit; (b) departure from 
the unstable equilibrium position; (c) system of primary and secondary dumb- 
bells to produce damping. 
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The net torque about the mass center C produced by gravity forces is* 
{F A — F B )a where the moment arm a is shown in Fig. 2(a). The gravity 
torque acts in such a direction as to diminish the angle 9. That is, it is 
a restoring torque which will rotate the dumbbell axis back to the local 
vertical. When the dumbbell becomes aligned with the local vertical, 
the moment arm a vanishes. Hence, the gravity torque becomes zero in 
the equilibrium position, = 0. However, due to the inertia of the masses, 
the dumbbell does not stop in its equilibrium position but continues to 
rotate past it, whereupon the gravity torque reverses its direction and 
acts to restore the dumbbell to the local vertical. This process produces 
an oscillation or "libration" about the local vertical which would con- 
tinue indefinitely if not damped out by some energy dissipating mech- 
anism. 

It is primarily the method of damping of the libration which dis- 
tinguishes the various gravity-gradient schemes from each other. One 
method of damping the librations requires the use of a second dumbbell. 
In order to understand the function of this second body we should point 
out that a dumbbell is also in equilibrium (that is, no gravity torque 
acts upon it) when its axis is perpendicular to the local vertical. How- 
ever, this equilibrium position is unstable in the sense that when the 
dumbbell deviates from the local horizontal by an arbitrarily small 
angle <p, the gravity torque (F A — F B )a acts in such a manner (see 
Fig. 2b) as to increase the angle <p, i.e., to drive the system away from 
its (unstable) horizontal equilibrium position. 

The inherent instability of a horizontal dumbbell may be used to 
design an efficient oscillation damper shown schematically in Fig. 2(c). 
In Fig. 2(c) the primary dumbbell AB is connected by means of a 
frictionless hinge to a secondary dumbbell A'B'. A spring is placed 
between the two dumbbells which keeps them crossed at right angles 
when the spring is not stressed. An energy dissipating device (repre- 
sented in Fig. 2(c) by a piston in a close-fitting cylinder) is placed 
between the two dumbbells so that any relative motion of the two bodies 
results in a loss of mechanical energy (mechanical energy converted 
into heat energy). When the main dumbbell is deflected through an 
angle 9 from the local vertical, it experiences a gravitational torque Ti 
which tends to restore it to the local vertical. At the same time, because 
of the spring, there is a tendency for the secondary dumbbell to be car- 
ried along through an angle 9' in the same direction as the angle 9, 
thereby producing a gravitational torque T 2 on it which tends to increase 

* Actually this is a slight oversimplification since F A and F B are not exactly 
parallel, but it adequately describes the main principle involved. 
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6' still further. The net effect is that the gravity torques T\ and T 2 tend 
to drive the two dumbbells in opposite direction, as shown in Fig. 2(c), 
thereby dissipating a relatively large amount of energy per cycle in the 
damping unit. The configuration shown in Fig. 2(c) will damp out oscilla- 
tions in the plane of the orbit; in order to damp out oscillations per- 
pendicular to the orbital plane, it is only necessary to add a second 
horizontal dumbbell which is perpendicular to the first one when its 
spring is unstrained. The two secondary dumbbells may be rigidly con- 
nected to each other and still provide damping in both planes of motion. 

2.2 Description 

The previous section describes the basic principles of the two-body 
system. In this section we discuss one method of reducing these principles 
to practice and describe the main features of all the major components 
of the system. These have reached a sufficiently high level of develop- 
ment for us to believe that they may be designed in detail for a specific 
experimental satellite. 

Fig. 3 shows schematically what an actual configuration might look 
like. The long vertical "mast" connects the satellite to an upper "deck 
assembly" which serves as the tip mass for the primary "dumbbell" 
and as the unstable body. The deck assembly consists of two crossed 
dumbbells which meet at a "hinge unit" that is connected to the mast. 
This "hinge unit" is actually a universal joint (or Hooke's joint) which 
also provides elastic restoring forces (springs) and energy dissipation 
devices (dampers). 

It should be mentioned that the deck assembly may be placed much 
lower on the mast rather than at its extremity as shown in Fig. 3, which 
illustrates a "high-deck" configuration. If the deck assembly is lowered 
to the vicinity of the satellite proper, the configuration will be referred 
to as a "low-deck" configuration (see Fig. 8 below for a schematic draw- 
ing of a low-deck configuration). From a dynamics point of view the 
two systems are identical. The high-deck configuration can be erected 
in a simple manner by the release of a single trigger which separates the 
satellite from the launching vehicle and simultaneously initiates exten- 
sion of all rods. The low-deck arrangement has the advantage of being 
much stiffer structurally (see Section III) than the high-deck configura- 
tion and is much less sensitive in its response to accidental misalign- 
ment of the various rods due to initial curvatures, thermal bending, 
or micrometeorite impacts which might cause plastic deformation. It 
has the disadvantage of a more complex erection sequence (to provide 
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Fig. 3 — A passive gravitational altitude control (PGAC) system configuration. 

clearance of the deck rods as they oscillate about the satellite body) 
and introduces the need to elevate a set of antennas above the height 
of the deck assembly in order to avoid electromagnetic difficulties. 



2.2.1 Extensible Rods 

A convenient way of erecting the rods in space is to use the STEM 
(Self-storing Tubular Extensible Member) units designed and developed 
by DeHavilland Aircraft of Canada, Ltd. These units consist of a beryl- 
lium copper tape (0.002 inch to 0.005 inch thick, and 2 inches to 5 
inches wide) which is stored on a drum prior to extension, in the same 
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manner as a carpenter's steel tape. However, unlike the carpenter's tape, 
the STEM tape has been preformed so that it tends to coil into a long 
straight tube when unwound from the storage drum, as shown in Fig. 4. 
The tape has a tendency to unwind spontaneously if not restrained from 
doing so. In fact, it is necessary to supply a governor which limits the 
extension speed to a safe level or else to provide a motor which drives 
the tape out at a controlled rate. Whichever is used, the motor or the 
governor mechanism, it could be located at the extremity of the deck 
rod to act as part of the necessary tip mass, as shown in Fig. 3. 

2.2.2 Damper Unit 

The damper unit will permit the deck a motion of two degrees of free- 
dom with respect to the mast in the manner of a universal joint. A pro- 
posed damper unit is shown in Fig. 5, where the two rotationally sym- 
metric housings arc rigidly fixed to one another with their axes crossed 
at 90°. The deck assembly is free to rotate about the axis of the upper 
housing while the mast is free to rotate about the axis of the lower hous- 
ing, thus providing the desired two degrees of freedom. To provide the 
required restoring torque, the deck assembly is fixed to a rotor whose axis 
is aligned with that of the housing by means of two fine torsion wires (or 
ribbons) as shown. These wires are maintained under suitable tension by 
means of the leaf springs at each end of the housing. This taut wire pro- 
vides the rotational restoring torque required and also serves to keep the 
rotor axis aligned with the axis of the housing. A slot is provided so that 
the connecting rod to the deck assembly may rotate through a total angle 
of 120° before bottoming on the end of the slot. Although it is not antici- 
pated that the rod will ever hit its stops except for rare periods of tum- 
bling (following injection, or collision with a micrometeoritc) one may 
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Fig. 4 — Extensible rod element. 
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TORSION WIRE 




Fig. 5 — Damper unit. 



design all stops so that they have no tendency to cold-weld in space. 
Similar stops are provided to prevent excessive lateral or axial motion 
of the rotor during periods of tumbling or during the launch phase. 

Damping is provided by means of one or more bar magnets fixed 
along a diameter of the rotor. These magnets have horseshoe-shaped 
pole pieces which enclose an annular disk of permeable material (e.g. 
cold-rolled steel) whose outer rim is fixed to the housing. A small gap 
always exists between the faces of the pole pieces and the permeable 
disk by virtue of the accurate elastic suspension, and even when the 
rotor is bottomed during launch or tumbling, the stops maintain a pre- 
determined clearance. As the rotor turns with respect to the housing 
because of satellite oscillations, the pole pieces rotate magnetic domains 
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in the permeable disk, thereby creating magnetic hysteresis losses. Mag- 
netic hysteresis losses are particularly desirable because they depend 
essentially upon the amplitude of oscillation rather than the frequency 
of oscillation and have been found to be effective at the very low libration 
rates, which are of the same order of magnitude as the orbital frequency 
(a six-hour period in the case of a 6000-nm altitude). 

Damper units of the type described above have been developed which 
provide the estimated damping torques required for 6000-nm orbits. 
They appear to be sufficiently rugged to withstand the launching en- 
vironment and weigh approximately two pounds for the complete damper 
unit. The feasibility of constructing dampers for higher and lower orbits 
has been demonstrated. Further details of the damper development 
program are described in Section VII. 

2.2.3 Packaging and Multiple Launching 

Design layouts indicate that there is no difficulty in packaging a 
stack of several satellites within the confines of a suitable rocket vehicle 
in such a manner that they will withstand rocket thrust and vibrations 
with a minimum amount of additional structure weight. The packaging 
arrangement is such that each individual satellite may be ejected with 
the required separation speed (see Section VI). The rod extension process 
may be triggered by the same explosive bolt mechanism which causes 
satellite ejection. One possible method of achieving this is indicated in 
Fig. 6 which shows how the rod constraints, which are needed during 
the launch phase, are automatically removed when the satellite is in- 
jected into orbit. 

2.2.4 Weight Breakdown 

For a 235-lb satellite body, operating at a height of 6000 nm, computer 
solutions based on the work of Ref. 4 indicate that a good design is 
achieved if the principal moments of inertia of body 1 (principal dumb- 
bell) are I\ = 3333, L> = 3333, h = 10 lb-ft-sec", and the principal 
moments of inertia* of body 2 (secondary rod configuration) are h = 
450, h = 1000, h = 1450 lb-ft-sec 2 . These moments of inertia may be 
achieved by using a 60-ft mast rod, four 40-ft deck rods, and the mass 
distribution shown in Table I. 



* The given values of 7 4 , I 5 and I e differ slightly from those given in the ex- 
ample in Ref. 4; this difference is the result of computer studies made after Ref. 
4 was submitted for publication. 
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Fig. — Packaging of satellites for multiple launching. 
III. ELASTIC VIBRATIONS VERSUS RIGID BODY MOTION 

The dynamics analyses 4 have been based on the assumption that the 
rods behave in an essentially rigid manner under the "gravity-free" 
conditions and extremely low librational angular speeds which prevail 
in the anticipated orbits; yet by ordinary earth-bound standards, the 
long thin rods would seem extremely flexible. To justify the assumption 



Table I — Mass Distribution 



Deck tip masses 
Deck rods 
Damper assembly 
Mast rod 
Mast motor 
Support structure 


29 lbs 
2\ " 

2 " 
4* " 
5* " 

3 " 


Total 


46 " 
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of rigid rods in the basic dynamics studies, it is necessary to estimate 
the bending and twisting deflections which occur in service. 

It is known that under normal circumstances the entire satellite will 
be oscillating, or librating, at certain well defined frequencies which 
are of the order of orbital frequency fi/27r. If any part of the system, 
which is supposed to behave in a rigid manner, happens to have a 
natural vibration frequency of the order of Q./2r, large deformations 
might occur, and the results of the rigid body dynamics analyses would 
be open to serious question. It is the purpose of this section to show that 
PGAC may be designed so that its smallest natural frequency is well 
above the orbital frequency 9./2v. 

It will be assumed in this section that the satellite has been aligned 
along the local vertical and that the elastic members of the system are 
undergoing small vibration with respect to a nonrotating* set of co- 
ordinate axes. In order to rigorously calculate the lowest natural fre- 
quency of such a system, it would be necessary to consider a set of 
coupled partial differential equations of considerable complexity. How- 
ever, it will be shown in Section 3.1 below that the mast rod may be 
considered to be perfectly rigid in the frequency range of interest. This 
result enables one to find the bending and twisting frequencies of the 
deck assembly in a relatively simple manner (involving ordinary rather 
than partial differential equations) as shown in Section 3.2. Finally, it 
is shown in Section 3.3 that the torsional mode of the mast rod has the 
lowest natural frequency of interest but that this frequency is still 
several times higher than the orbital frequency. 

3.1 Rigidity of Mast Hod and Influence of Distributed Mass oil Natural 
Frequencies 

Let us consider the bending vibrations of a long beam, of length L, 
flexural rigidity EI, and mass p per unit length, whose ends carry two 
relatively large tip masses, Mi and Mi , but which are otherwise un- 
constrained, as shown in Fig. 7. Since the tip masses are so great com- 
pared to the beam mass il/ 3 , the tips can never move too far from their 
equilibrium position (in comparison with the midpoint of the beam). 
Therefore, the principal mode shape must look somewhat as shown in 
Fig. 7 with nodal points very near the ends. We can thus consider the 
problem equivalent to that of a beam, of length L' ~ L, whose ends are 
fixed against displacement and more or less fixed against rotation, de- 
pending upon the constraints provided by the tip masses. In any case, 

* This is equivalent to neglecting (ho curvature of the orbital path. 
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Fig. 7 — Mode shape for light, beam carrying heavy tip masses. 

the circular frequency p of the fundamental mode is given (Ref. 8, 
pp. 324-339) by 

(1) 



v = (C/L')WEI/p 

where C is a numerical coefficient which depends upon the degree of 
constraint at the nodes. 

Now consider two limiting cases of constraint between the mast 
and deck assembly. If the springs which connect the deck and mast are 
extremely soft, no appreciable bending moment can be transmitted to 
the mast from the deck, so the connection point may be treated as a 
hinged or simply supported end. If the tip mass at the other end has 
negligible moment of inertia, that end may also be considered as simply 
supported and the coefficient 8 C = 3.14 (hinged-hinged beam); but if 
the tip mass has an appreciable moment of inertia, the end may be 
considered clamped, in which case C = 3.93 (hinged-clamped beam). 
Another limiting case arises if the connecting spring is extremely stiff, 
in which case the deck assembly with its very large moment of inertia 
is almost rigidly fixed to the mast rod and essentially prevents rotation 
of the connected end of the mast. Tn this case C = 3.93 (clamped-hinged 
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beam) or C = 4.73 (clamped-clamped beam) accordingly as the tip 
mass has negligible moment of inertia or infinite moment of inertia. 

We thus see that C can only undergo moderate variations despite 
extreme changes in the manner of end support. If one adopts the most 
conservative point of view and considers C = 3.14, (1) predicts that 
for a typical mast rod (L « U = (500 in, EI = 38,300 lb in 2 , p = 
1.51 X 10 -5 lb sec" in - ") the fundamental mode of vibration has a fre- 
quency of p = 1.37 rad/sec, which is certainly well above the orbital 
frequency Q = 0.273 X 10 - ' rad/sec for a 6000-nm orbit. Thus, we see 
that the influence of the distributed rod mass cannot play an important 
role in vibrational motions at the low frequencies of interest, and the 
mast rod may be considered rigid. 

3.2 Natural Frequencies of Deck Assembly 

Having shown that the mast rod is practically rigid, one may con- 
sider the satellite and mast rod a single rigid body upon which is mounted 
the deck assembly via the flexible joints of the damper unit. Because 
the deck tip masses are so great compared to the mass of the deck rods, 
one may neglect entirely the deck rod mass and treat the problem 
according to the standard "lumped mass" point of view. In particular, 
if one restricts attention for the time being to the low-deck configuration 
illustrated in Fig. 8, and notes that the center of mass* of the entire 
system lies fairly close to the plane of the deck assembly, one may 
introduce a further simplification by considering the deck assembly to 
be oscillating about a fixed point where the two hinge axes are assumed 
to cross. 

In setting up the equations of motion, we shall use D'Alembert's 
principle, wherein each moving mass Mi is thought of as loading the 
structure by a system of "inertia forces": Xi = —Mj&i , F, = — Mm , 
Zi = — Mi&i , parallel respectively to the axes of x, y and z; and the 
rods are loaded by "inertia torques", 7\ = —Infii . Superscript dots 
denote differentiations with respect to time t, in the usual Newtonian 
notation. To provide a more flexible and symmetrical notation, we shall 
frequently speak of the generalized displacements, </, , and generalized 
forces, Qi , which are related to previously defined quantities as in 
Table II. In this tabulation, generalized masses ?«, have been defined for 
future reference. 

Each independent deflection </, can be found as a function of the 



* The center of mass lias been assumed to be unaccelerated in inertia! space and 
for our purposes may be considered fixed. 
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Fig. 8 — Low-deck configuration. 
inertia forces applied to the system 



(2) 



The terms dgt/dQj are called influence coefficients and are functions of 
the elastic constants and dimensions of the system. Each influence 

Table II — Generalized Displacements, Forces and Masses 
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Table III — Low-Deck Configuration 



Frequency />„ (rail/sec) 



p,= 0.597 X 10"* 
p« = 0.G43 X 1<T 3 



p 3 = Pi = Pb = 0.401 X 
]H = pi = V* = 0-SOO X 
p« — P>u= 0. 1373 
pn = p\< = 0.2025 



io-> 

10' 



Most Significant Vibration 



Oscillation of deck assembly about x (roll) 

axis 
Oscillation of deck assembly about y 

(pitch) axis 
Bending of deck rod parallel to X-axis 
Bending of deck rod parallel to y-axis 
Twisting of deck rod parallel to x-axis 
Twisting of deck rod parallel to ?/-axis 



coefficient may be found by an elastic analysis of a statically determinate 
structure; for the sake of brevity, these coefficients will not be explicitly 
written out here. 

Equation ( 1 ) may be written in the form 



£[&*-«* 



= 



(3) 



where 5„ = 1 for i = j, and 8,,- = for i ^ j. The symbol Qj represents 
the generalized inertia force, — wi,-g> , and m, represents the generalized 
mass (or moment of inertia) defined in Table II. The equations of 
motion are thus given by 



A 17*. 

h L W 



(—mjqj) — Sijqi 



= 0. 



(4) 



In order to solve the system of differential equations, we may assume 
that 

q t = Ai cos (pt — ip) (5) 

where Ai , p and <p are as yet unknown quantities. If one substitutes 
(5) into (4) and follows the standard procedure, one finds a twelfth- 
degree equation in 1/p 2 with twelve solutions, l/p„ 2 for n = 1,2, • • • , 12. 
For these twelve (not necessarily distinct) values of p, (4) is satisfied 
and the assumption of (5) is justified. Equation (5) shows that the 
terms p n represent the circular frequencies of the so-called natural 
modes of vibration. Upon the introduction of suitable numerical values, 
one finds the twelve natural frequencies pi • • • p& , listed in Table III, 
for a typical configuration designed to orbit at 0000 nm. 

It may be seen from Table III that the lowest natural frequencies 
arc those corresponding to the oscillations about the hinge-spring axes.* 

* The frequency of the "rigid body" oscillations of the deck assembly about the 
pitch and roll hinge axes have been made intentionally close to the libration fre- 
quency in order to provide good damping. All other natural frequencies must be 
kept well above these values to avoid undesired resonances. 



2212 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1903 

These frequencies differ from those which would be obtained with per- 
fectly stiff deck rods by one part in 5000. The bending of deck rods has 
the next highest natural frequencies, which are about 70 times the pitch 
spring frequency. This indicates that excitation at a libration frequency 
(%0.5 X 10 -3 rad/sec) would not cause very large unwanted deflections 
any place in the structure, and that the assumption of rigid rods in the 
dynamics analysis is well justified for the low-deck configuration. 

Although a complete vibration analysis of the high-deck configuration 
(shown schematically in Fig. 3) has not been made, there is no reason 
to believe that the natural frequencies of vibrational modes dominated 
by bending action will differ by orders of magnitude from similar modes 
in the low-deck configuration. 

3.3 T(rrsional Oscillations of Mast 

On the other hand, it is to be expected that the frequency of torsional 
vibration about the mast axis will be considerably less for the high-deck 
configuration. As a first approximation, one may neglect the bending 
deformations of the deck rods and consider the system shown in Fig. 3 
as a long rod of torsional constant K m (K m = torque per unit twist angle) 
separating two rigid bodies whose moments of inertia are h and la , 
respectively. The angular frequency of natural oscillation is given (Ref. 
8, p. 12) for such a system by 



V = 



^H!I 



where the approximation follows from the fact that Id^> /& . For a typical 
case of interest, one would find a torsional oscillation frequency, for 
the high-deck configuration, of the order of 

p = 0.0048 rad/sec. (7) 

This value should be compared with a libration frequency (in a so- 
called higher roll-yaw mode) of p Ty ^ 0.00049 rad/sec. If a somewhat 
more refined analysis is made, which takes into account the elasticity of 
the deck rods, the improved value of p differs insignificantly from the 
value given by (7) . Although this value is smaller by an order of magni- 
tude than the corresponding frequency of the low-deck configuration, it 
is still about ten times greater than the largest libration frequency. Some 
other comparisons between high- and low-deck arrangements have 
already been discussed in the introduction to Section II. 
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IV. STRESS AND DEFLECTION ANALYSIS OF RODS DURING EXTENSION PHASE 

Since a satellite cannot be injected into orbit with absolutely zero 
angular velocity, in inertial space the tip masses on the extending rods 
will tend to cause bending and twisting of the rods during the process 
of extension. No attempt will be made to examine this problem in full 
generality, but two important representative cases will be considered. 
In both cases only the high-deck configuration is considered, since the 
low-deck configuration would seem to be at least as strong as the high- 
deck configuration. 

Experiments 910 have demonstrated that a properly designed spring 
arrangement is capable of injecting satellites into orbit with tumbling 
rates below 1 rpm prior to rod extension. It is shown in this section that 
such rates do not cause excessive stresses or deformation in the rods 
during the extension process. 

4.1 Tumbling 

The satellite is idealized as shown in Fig. 9 and is assumed to be 
tumbling at time I = with angular speed w about the body axis x, 



-DECK ROD 




CM. OF SATELLITE 



M b ,Ib 



SATELLITE 
BODY 



Fig. 9 — Schematic diagram of satellite during extension. 
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which is parallel to the deck rods carrying tip mass Mi . At time t = 
both mast and deck rods begin to extend with speed v, as indicated in 
the sketch where r represents the initial distance between the deck tip- 
masses and the mast axis, and R equals the initial distance, measured 
along the mast axis, between the satellite body and the deck assembly. 
All mass in the deck assembly exclusive of the rods and tip masses is 
considered to be concentrated at the tip of the mast (point A) in the 
rigid body labelled "deck structure" in Fig. 9. The mass of the "deck 
structure" is denoted by Ma , and its moment of inertia about a ccntroidal 
axis parallel to x is denoted by I d ; similar expressions for the satellite 
body are denoted by Mb and h , respectively. 

From Fig. 9 it is seen that the instantaneous distances r, Ri , and Ro 
are given by 

r = r + vt 

(2Mi + 2M 2 + M d )(R + vt) 



Ri = 



R-2 = 



2(M 1 + M 2 ) + M d + M b (8) 

M b (R + vt) 



2(Mi + Mi) + M d + M b 

and the instantaneous moment of inertia / of the entire system about 
the rc-axis (passing through the instantaneous center of mass) can be 
shown to be 

I(t) - I b + la + M(R + rf) 1 + 2M,(r + vt) 2 (9) 

where M is defined by 

, 7f = M b (2Mi + 2M 2 + Mg) . . 

2Mi + 2M 2 + M d + M b ' V ; 

The initial value of / is denoted by 7 and is found from (9) by setting 
t = 0. 

We shall now assume that: (i) the mass of the rods is negligible; 
(m) the hinge connection between the mast and deck assembly is rigid; 
(in) the rods do not bend or twist (until further notice); and (iv) the 
mass center of the system is moving through inertial space with constant 
velocity. Under these assumptions one may apply the principle of con- 
servation of angular momentum* to find the angular velocity to and 
acceleration a> in the form 



* Although unguliir momentum is not strictly conserved in the presence of 
gravity torque, it can be shown that this effect is not significant. 
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Wo/(] 



w = 



— coq/q dl 

P dt' 



(11) 



Equations (11) together with (9) fully specify the angular velocity and 
acceleration during the entire extension phase. The absolute acceleration 
a'' of any point P in the system may be found from the vector equation 



a" = 5 J? + * x p + 2o> X £ + o> X (» X p) 



8t 2 81 



81 



(12) 



where p is the position vector of the point P measured from the origin 
shown in Fig. 9, or in terms of the unit vectors i, j, k along the body 
axes: 



By definition: 



p - pA + p v j + p,k. 

Sp/5/ = pA + p,j + 0«k 
5 2 p/5< 2 = pj + p„j + p\k 
to = o>i; 5o>/8/ = wi. 



(13) 



(14) 



Applying (12) successively for the five points A, B, C, D, and E shown 
in Fig. 9 one may find the components of acceleration d? = a x i + 
a„ P j + a 2 P k indicated in Table IV. 

Table IV, together with (11), gives the absolute acceleration of all 
the tip masses. The D'Alembert forces acting on masses at points A, B, 
C, D, and E are respectively -Mj& A , -M&", -MiSi c , -il/ 2 a D , -il/ 2 a*. 
The bending moment M mx at any point z along the mast is given at any 
time by 

M mt = -[M d a y A + Miiaf + O + MiW + a u E )](R 2 - z) 

+ M 2 r(a 2 E - a 2 ") - I d u. (15) 







Table IV 


— Components of Acceleration 


p 

(Point) 


Px 


Pa 


Pu 


«,' 


««* 


a z p 


A 
B 
C 
D 

E 



r 

— r 








r 


R, 

Ri 
Ri 
Ri 
R, 









-(2fl2W + RiOl) 

- (2fijw + R«.io) 
-(2/^*0 + R-iu) 
rur — 2R&> — rt»w 
-(no 2 + 2R-& + R t u) 


-R#>* 
-flaw 1 
-R#fi 

-(flaw 2 + 2rw + rc6) 
-R#) 2 + 2fu + rd 
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Table V — Parameters 



M b = 286/32 lb sec 2 /ft 


h = 10 or 20 lb sec 2 ft 


Mi = 8/32 lb sec»/ ft 


Id = 0.0845 lb sec 2 ft 


Mi = 4/32 lb sec 2 /ft 


Ro = 2 ft 


Mi = 9/32 lb sec 2 ft 


r = 0.75 ft 


EI = 13.4 lb ft 2 (deck rod) 


EI = 297 lb ft 2 (mast rod) 



Similar expressions are readily written down for the bending moments 
at various points along the deck rods but will be omitted here for the 
sake of brevity. Deflections are found by noting that a tip force F 
produces a lateral tip deflection A F = FL /(ZEI) for a cantilever of 
length L and bending stiffness EI. Similarly, a tip couple M produces 
a lateral tip deflection of amount A. w = ML 2 /(2EI). The net deflection 
is found by superposition. For a GOOO-nm satellite similar to the one 
described in Section II, the parameters shown in Table V were used. 
An investigation of the complete extension history shows that for an 
initial tumbling rate of 0.1 rad/sec (^1 rpm) and an extension rate of 
v = | ft/sec, the maximum stresses occur early in the process and decay 
rapidly thereafter; i.e., maximum moments occur before the rods have 
extended a distance of 2 ft. The maximum bending moments (which 
occur at the cantilever root) and the corresponding tip deflections 
(expressed as a fraction of rod length at the instant of maximum loading) 
are given in Table VI. Published data 11 indicate that short lengths of 
the mast rod could sustain a bending moment about a hundred times 
greater than the maximum value indicated in Table VI, and the deck 
rods could sustain a value about 30 times larger than the greatest 
tabulated value. Thus, there appears to be no "stress" problem due to 
tumbling. 

4.2 Spinning 

Assumptions (i) to (iv) of the previous section will be retained. 

If the entire satellite spins about the mast axis with angular speed 



Table VI — Bending Moments and Deflections for v = \ ft/sec; 



o>o = 0.1 rad/sec 



Satellite Moment 
of Inertia 


Maximum Bending 
Moment, Mast 


Maximum Bending 
Moment, Deck 


Mast Tip 
Deflection 


Deck Rod Tip 
Deflection 


0b sec= ft) 
10 

20 


(ft-lb) 

0.072 
0.113 


(ft-lb) 

0.018 
0.029 


(Fraction of rod 
length) 

0.0003 
0.0006 


(Fraction of rod 
length) 

0.0012 
0.0025 
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wo at time / = 0, the deck assembly will acquire an angular speed B d 
and the satellite body will rotate at speed B b . Conservation of angular 
momentum* requires that 

[77 + 2(M, + J/--)(r + vt)% + I A 

= [// + 2(Mi + M*)ro + h]m (16) 

where h' denotes the moment of inertia of the "deck structure" about 
the mast axis. The torque required to produce a unit relative angular 
displacement between the satellite body and the deck is proportional 
to (R — RoV 1 and may be represented in the form k/(vt), where k is 
the torsional rigidity of the mast for unit length. In other terms, the 
torque on the mast at any time is given by (0 d - 0,,)k/vt; this torque is 
applied directly to the satellite body, so one may write 

IA = k(6 d - 6 b )/vt. (17) 

Equations (16) and (17) represent a third-order system of linear dif- 
ferential equations with time-dependent coefficients and with initial 
conditions of the form d (O) = 0; 6 (O) = 0; 6 (O) = w . These equations 
have been integrated numerically to provide the complete response of 
the system during the extension phase for various sets of parameters. 
The solutions indicate that torsional stresses do not become excessive 
at any time, although the satellite body might rotate, relative to the 
deck assembly, by as much as 10 revolutions if an extension speed of 
v = \ ft/sec is used, and o> is as high as 2 rpm. The bending stresses 
and deflections produced in the deck rods, under these conditions, are of 
the same order of magnitude (very safe) as found in the foregoing 
section on "tumbling." With a nonspinning final-stage vehicle it is 
unlikely that the initial spin rate co will reach a value as high as 1 rpm. ' 

4.3 Umbrella Effect 

If all rods are being extended simultaneously, for the high-deck 
configuration the tip masses on the deck rods will continue to move 
parallel to the mast axis at the termination of the mast extension 
phase. This motion will continue until the cantilever bending of the 
deck rods has converted the tip mass kinetic energy into stored elastic 
energy. This effect will be referred to as the "umbrella" effect. To com- 
pute the maximum tip deflection 5,„ and the maximum root bending 
moment M ,„ in the deck rods, it should be observed that a lateral tip 



* Effect, of gravity torque is neglected here, as in Section 4.1. 
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force P will produce a tip deflection 8 m = PL Z /SEI, where L is the beam 
length and EI is the fiexural rigidity. The stored strain energy U is 
equal to (P/2)8 m , which may be expressed as U = SEI8 m '/2L' by 
virtue of the linear relationship between P and 8. When the stored 
energy U is equated to the initial kinetic energy (Mv /2) of a tip mass 
M which moves at speed v, one finds the tip deflection 

8 m = v[ML 3 /3EI)}\ (18) 

The root bending moment M m is found by multiplying the load P = 
(BEIdm/L) by the beam length L to give 

M m = v[SEIM/L\\ (19) 

For a typical deck unit with the following parameters, EI = 13.4 lb 
ft 2 , L = 50 ft, M = (10/32.2) lb sec 2 /ft, one finds that for sufficiently 
small values of v 

8 m (ft) = 0.983 v (ft/sec) 

M m (ft-lb) = 0.500 v (ft/sec). 

Thus an extension speed of v = \ ft/sec would produce a bending 
moment of about ;} ft-lb, which is less than the approximate allowable 
value of 1 ft-lb. The corresponding tip deflection of about \ ft is suffi- 
ciently small so that the linear bending theory used is adequate. If one 
wishes to find the maximum extension speed which produces a root 
bending moment below 1 ft-lb, it is necessary to consider a nonlinear 
beam theory which allows for large slopes in the deflected beam shape. 
An approximate treatment of this problem indicates that an extension 
speed of about 1.5 ft/sec would result in a root bending moment of about 
1 ft-lb. Therefore, if one wishes not to exceed the load-carrying capacity 
of the deck rods, it is essential to keep the extension speed well below 
1.5 ft/sec, or else to extend the deck rods after the mast has been fully 
extended. 

Similar considerations show that reasonable differences in the exten- 
sion speeds of the various deck rods result in tolerable loads on the mast, 
for practical configurations. 

It should be noted that when the oscillating deck masses slam down- 
ward, they load the mast axially and tend to produce Euler-type buck- 
ling. It may readily be shown that, so long as the mast extension veloc- 
ities are kept small enough to prevent overloading of the deck rods, the 
mast will not buckle for the configurations of interest. 
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V. THERMAL LOADING 

The rods used in the proposed design consist of long split overlapping 
tubes which will experience temperature gradients due to solar heating. 
These temperature gradients will cause the rods to bend in such a way 
that the illuminated side becomes convex. 

In this section, the lateral deflections of the rods due to solar heating 
will be calculated. In Section VIII it will be shown that these deflections 
have a minor influence upon the pointing accuracy of the PGAC System. 

5.1 Temperature Distribution 

It will be assumed that the rod is sufficiently long so that end effects 
may be ignored ; hence the temperature distribution will not vary with 
length along the rod. Since the heat input depends upon the angle be- 
tween the collimated solar rays and the axis of the rod, it is implicit in 
the above statement that the thermally induced curvature of the rod 
axis is small; this necessary requirement will be verified a posteriori in 
a numerical example. Confining attention to a unit length of rod as 
shown in Fig. 10, it may be verified that the cosine of the angle between 
the solar rays and the normal to a surface element located at an angle 
0, measured from the outer edge of the tape, is given by 




° N,T Let 



Fig. 10— Unit length «>f split overlapping tube illuminated by the sun. 
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sn = sinp, cos (0„ — 0) (20) 

where s is a unit vector pointing to the sun, n is a unit surface normal, 
ip, is the angle between the solar rays and the tube axis, and 6 a is the 
angular distance from the outer free edge of the tube to the plane formed 
by s and the tube axis. The heat input per unit time on a unit area of tube 
surface is given by 

g s = a S sin <p a cos + (0 — S ) = a S cos + (0 — 0„) (21) 

where 

a = absorptivity for solar radiation 

S = solar constant (442 Btu/hr ft 2 ) 

S = S sin <p s = effective solar constant 

cos + x = \ (cos x + I cos x |) (half -rectified cosine wave). 

In general, a small element of the tube of arc length rdd (where r is the 
tube radius) gains heat q 3 rdd, in unit time, due to solar heating; the 
element also gains heat q e rd& by conduction and loses heat q e rdd by 
emission of radiation. It will be assumed that the overlapping layers do 
not have an appreciable area in mutual contact. This idealization is 
useful because of the random nature of the actual contact areas and the 
uncertainties in the contact pressure and in the associated surface heat 
transfer coefficients. Any heat conduction which does occur between 
overlapping layers will tend to reduce temperature gradients and 
alleviate the thermal bending effect; thus, the neglect of such effects 
leads to a conservative analysis. Since the walls are very thin, it is per- 
missible to assume that the temperature varies only in the circum- 
ferential direction, so that Fourier's law leads to the result 

q c = (Kh)(d 2 T/rW) (22) 

where k = thermal conductivity, h = wall thickness, and T = absolute 
temperature. The heat loss by radiation is given by the Stefan-Boltzmann 
law: q e = toT , where e = hemispheric emissivity at temperature T, 
a = Stefan-Boltzmann constant (1714 X 10"' 2 Btu/hr ft 2 (°R) 4 ). For 
simplicity, the effects of internal radiation will be neglected, so that the 
above expression for q c is valid only for 5^ ^ 2t. The inclusion of 
internal radiation effects would result in reduced temperature gradients, 
thereby reducing the thermal bending; thus, this assumption is also 
conservative. 

Upon summing up the three contributions to the thermal balance, 
one finds that 
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g - (S*)l* = -(^ cos + (0 - Oi £ £ 2* (23) 



.i2rn 



g! = 0; > 2x. (24) 

Because no appreciable amount of heat can be radiated over the 
narrow faces (of area h per unit length) at the edges where = and 
= 0,„„x , one may write the boundary conditions in the form dT/dd = 
at = and = max , and observe that both T and dT/dd must be 
continuous at = 2?r. Equation (24) implies a linear temperature 
distribution in the range > 2tt, but since dT/dd vanishes at the edge 
= n .nx , the temperature must be constant in the range > 2ir. In 
addition, continuity of dT/dd requires that dT/dd = at = 2ir. 
Thus, the temperature distribution may be found by solving the non- 
linear differential equation (23), subject to the boundary condition 
dT/dd = at = 0, and at = 2ir. This problem, except for the boundary 
conditions, is similar to the problem treated by Charnes and Raynor" 
of a continuous (nonsplit) tube. Following their treatment, one may 
linearize (23) by writing 

T = r„+ r(d) (r«7' ) (25) 

where T is the mean radiant temperature defined by 

7'o = [aS/irw]*. (26) 

Upon substitution of (25) into (23), one finds 

d\/dd- - v'r = -^ cos + (0 - 0.) + X7o 4 

dr/dd = at = 0, and at = 2tt 

where 



rj 



(27) 



p = 2vWj x = (<r*r 2 /*h); /3 = aSr 2 /Kh (28) 

Equation (27) may be solved by a number of standard procedures (e.g., 
use of Duhamel integral or of Laplace transform) which will be omitted 
for the sake of brevity. For a typical beryllium-copper rod, the pertinent 
dimensions are 

r = 0.225 in, h = 0.002 in, k = 65 Btu/hr ft, d max = 3tt. 

Representative values of absorptivity and emissivity, calculated by 
integration of monochromatic reflectivity measurements, are: 

a = 0.8, e = 0.3 
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Fig. 11 — Temperature distribution in split overlapping tube (neglecting 
internal radiation and radial conduction effects). 

where e corresponds to a temperature of 7'o = 684°R. If the solar rays 
are truly normal to the rod axis, one may use S = S = 442 Btu/ft 2 hr 
and if the rod is oriented with 6 S = 7r/2, as shown in Fig. 11, the tem- 
perature distribution will be as shown in the figure, where t = T — T 
is plotted radially outward from the tube surface for positive values and 
inward for negative values of t. For this example, the temperature drops 
continuously from T = 700°R (240°K) at 9 = to T = (i07°R (207°F) 
at = 3tt. 

5.2 Thermal Bending 

Following the method used by Timoshenko and Goodier 13 for bending 
of a beam of rectangular cross section, one may show that the curvatures 
developed in the z-x plane and the z-y plane are given, respectively, by 

K X = My/Ely, K U =MJEI X (29) 

where E is Young's modulus, 



2 dA 



h = f ifdA, I y = f x 2 
\I Z = Ea I rydA, My = Ea I rxdA. 



(30) 
(31) 
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In the above expressions, a represents the coefficient of thermal 
expansion (a = 9.4 X 1(T 6 °R~' for numerical example); x and y are 
measured from the centroid of the cross section which is located a 
distance e from the axis of the tube, as shown in Fig. 1 1 ; x and y have 
their origin at the centroid; dA represents an element of area; and the 
integration is made over the entire cross section. 

With the temperature distribution shown in Fig. 11 (corresponding 
to the numerical data given above), one finds that the curvature k x is 
negligible, but k u is found to be 

Ky = 8.18 X 10 3 ft -1 (R u = 1/k v = 314 ft). 

It may readily be shown that one end of a rod of length L bends 
through an angle A^ = kL with respect to the other end and deflects 
through a lateral distance of 8 = (%)L~/R. The maximum angular and 
lateral deviations for a 50-ft length of rod are thus seen to be A^ = 9.1°, 
5 = 4.1 ft. A similar calculation shows that 5 = .3.5 ft for a rod of r = 
0.45 in and h = 0.005 in. 

In view of the conservative nature of the heat transfer analysis used, 
one may estimate that the actual values of slope and deflection could 
easily be less than half of the computed values. In any case, the deflec- 
tions do not cause excessive misalignment from the local vertical (see 
Section VIII), and the slopes are sufficiently small to justify the initial 
assumption that the heat input and temperature distribution do not 
vary appreciably along the axis of the rod. 

VI. SPRING DESIGN FOK MULTIPLE LAUNCH 

When several satellites are launched from the same rocket vehicle, 
it is necessary that they be injected with different velocity components 
along the orbit trajectory; otherwise, all the satellites will have the 
same period and will appear to be "bunched" together when viewed 
from the ground. The velocity increments required to "minimize" the 
undesirable effects of bunching are discussed in Ref. 14, where it is 
indicated that a relative speed of about 12 ft/sec between the slowest 
and fastest satellites is desirable for the case of four satellites in a single 
orbit at 6000 run. Similar conclusions were reached in unpublished work 
at Bell Telephone Laboratories for the case of three simultaneously 
launched satellites. 

Velocity increments of 12 ft/sec are readily achieved by mechanical 
springs. In this section we shall consider the use of ordinary helical 
springs and of the so-called conical disk-spring (sometimes called a 
Belleville spring) shown schematically in Fig. 12. 
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Fig. 12 — Belleville spring showing stress distribution. 

It is easily shown 15 that the strain energy per unit volume Uk stored 
in a close-coiled helical spring with a narrow circular cross section can be 
expressed in the form 

u h = {\)tJ/G (32) 

where t„, is the maximum shear stress in the spring and G is the shear 
modulus. 

In the case of a Belleville spring, the stresses are distributed 16 in an 
approximately linear manner over the cross section, as shown in Fig. 12, 
if the inequality (b — a) /a « 1 is satisfied and only small deflections 
are permitted. Under these conditions, it is easy to show that the strain 
energy per unit volume n b is given by 



n, = (IW/E 



(33) 



where a m is the maximum tensile stress in the spring and E is the modulus 
of elasticity. The relative energy-storing efficiencies of helical and Belle- 
ville springs may be found from (32) and (33) in the form 

Wb = (2\(G\(a» 
u h \ZJ\e)\t„ 



2* 



(34) 



where use has been made of the well-known relationship (Ref. 15, p. 60) 
(E/G) = 2(1 -f- v) and of the fact that Poisson's ratio v is very close 
to ^ for most structural metals. 

If the spring material is to be used most effectively, the stresses r m 
and ffm should be practically equal to their respective values at the elastic 
limit. It is seen from (34) that the relative efficiency of Belleville springs 
versus helical springs depends upon the ratio of (<r m /T m ) at the elastic 
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limit. This ratio depends upon the criterion of elastic failure which 
governs the material. For example, a relatively ductile metal tends to 
yield (Ref. 1G, Section 82) when either the shear stress or the octa- 
hedral shear stress reaches a critical value ; for such materials it may be 
shown that <r m « 2r m . For materials with little ductility, failure gen- 
erally occurs by fracture and <r m « r m (Ref. 17). Therefore, (34) pre- 
dicts that if the material is stressed up to its useful limit: 

— ^ 1 , for ductile materials 
u h 

— « - , for brittle* materials. 
Mft 4 

In other terms, both Belleville and helical springs require the same 
volume of any given ductile metal to store equal amounts of energy; 
but a Belleville spring can store only £ the energy stored in an equal- 
volume helical spring made of the same relatively brittle material. 
Although there is a tendency towards weight saving in the use of a 
helical spring made of a relatively brittle material, it may well be that 
practical geometric considerations, reliability, and the reserve strength 
of ductile metals would lead one to the choice of a Belleville spring. 

To show that reasonable spring weights are required for the present 
application, let us equate the strain energy in the spring to the kinetic 
energy (?)(W sat /g)v 1 required to impart a separation speed v to a 
satellite of weight W au t(g = 386 in/sec 2 ). If the volume of spring material 
is denoted by 7 S „ , (33) leads to the result 



U b - u b \ BP - _ _ (35) 

If one uses the relationship W BP = wV BP , where W BP is the total weight 
of the spring and w its specific weight, (35) leads to the conclusion that 

W Bp = ir. ttt — I— 5 • (3(3) 

Q v»7 

Equation (36) is a compact expression for the weight of a well designed 
Belleville spring or of a helical spring (for a material with o m W 2r m ). 
The material influences the spring weight only through the ratio 
(aJ/wE), which may be interpreted as twice the elastic energy stored 
in a unit volume of the material when uniformly stressed at its maximum 

* The words ductile and brittle are used in the sense that there either is or is 
not an appreciable amount of plastic flow between yield and fracture. 
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Table VII — Material Properties 



Material 


E 


ai 


*M 


■wE/Om- 


(i) 4130 steel 

(HT 200,000 psi) 
(u) Aluminum alloy 

7075T(G) 
(Hi) Titanium Alloy 
(Ti-0Al-0V-2Sn) 


(psi) 
29 X 10 s 

10.4 X 10" 

17 X 10" 


(lb/in») 
0.282 

0.101 

0.162 


(psi) 
175,000 

73,000 

190,000 


(in- 1 ) 
2.67 X 10- 4 

1.97 X 10-* 

0.763 X 10-« 



allowable value of «r* . Table VII shows (wE/aJ) for some typical 
materials of interest. Thus, if one wished to impart a velocity of v = 12 
ft/sec to a satellite weighing W sat = 280 lb, (36) shows that with the 
three materials described above, the spring weight W BJ > would be 12.1 
lb, 8.8 lb, and 3.4 lb for materials (»), (») and (Hi), respectively. 

VII. damper unit 

The damper unit was described in qualitative terms in Section 2.2.2. 
In this section it will be shown in what respects the PGAC damper 
differs from other dampers that have been proposed in the literature, 
and how the damping torques and spring torques must be chosen in 
order to meet the system requirements outlined in the Introduction, 
Section I. The hardware development program for the damper units is 

also described. 

Since the means of damping libration motions is perhaps the single 
most important feature which distinguishes the various attitude control 
systems that have been proposed by several authors, it would seem 
worthwhile to indicate the various methods that have been considered. 
These fall under two main categories: (a) velocity-dependent damping, 
and (b) amplitude-dependent damping. 

In the first category, one finds schemes which depend upon viscous 
fluids 2 ' 51819 or eddy currents. It has not been demonstrated that practi- 
cal difficulties concerning seals, viscosity, temperature and adverse 
rheological effects have been overcome in lightweight systems utilizing 
fluids. Calculations have shown that effective eddy-current damping 
requires a considerably greater weight of material than does the magnetic 
hysteresis unit under discussion. 

In the second category of damping methods, the energy loss per cycle 
is independent of velocity but depends only upon the amplitude of mo- 
tion. Included in this category are methods based upon Coulomb fric- 
tion, internal friction, 36 and magnetic hysteresis, as described in Section 
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II of this paper. Coulomb friction (also called dry sliding friction) de- 
pends upon physical and chemical surface properties which are notori- 
ously hard to control under conditions of high vacuum and thermal 
cycling; it is also difficult to control the normal force between the sliding 
bodies, which greatly influences the level of friction. Solid internal fric- 
tion, which depends upon energy losses developed in the microstructure 
of the material, is quite temperature-dependent but does not depend 
upon unreliable surface properties and should not be unduly influenced 
by high vacuum. 

In addition to the virtues of velocity independence, insensitivity to 
surface conditions, and lack of rubbing parts, the magnetic hysteresis 
damper proposed here has been shown to exhibit relative insensitivity 
to wide temperature fluctuations. 

7.1 Spring Constants 

It is shown that because the deck oscillates about an unstable position 
of equilibrium, it is necessary to satisfy certain stability criteria. This 
requires that the torsional spring constants fci and ko exceed certain 
critical values fci* and !;■>* given by 

kf/htf = 0.625 (roll) 

(.37) 
h?/U? = 1.3 (pitch) 

for the satellite specified in Section 2.2.4, or, for 6000-nm altitude, 

ky* = 0.155 X 10~ 3 ft-lb/rad 

/,-,* = 0.324 X 10~ 3 ft-lb/rad. 

It is also found that k cannot be too large, since the two-body system 
becomes so stiff at large k that very small relative displacements between 
the two bodies are developed and the energy dissipation due to amplitude- 
dependent damping is reduced. To guarantee stability, it has been 
decided to keep Ai and A - 2 at least 10 per cent above their critical values. 
Computer studies indicate that the variation in damping time is rela- 
tively small in the range: 

A-! = 0.175 X 10~ 3 to 0.375 X 10~ 3 ft-lb/rad, 

{h/htf = 0.7 to 1.5) 

3 •, (38) 

fc 2 = 0.36 X 10 to 0.76 X 10~ 3 ft-lb/rad, 

(h/I-S? = 1.45 to 3.06). 
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From computer solutions it has been noted that if k is smaller than k*, 
both the satellite body and the deck body will oscillate about cocked 
equilibrium positions. This verifies the stability criteria for the spring 
constants. If the spring constants are much larger than the maximum 
values given in the above ranges, the relative angular displacements 
become very small and little energy dissipation occurs. 

7.2 Damping Torque 

Damping torque is produced by rotational hysteresis losses obtained 
from relative displacement between a magnet, fixed along a diameter 
of the rotor, and an annular thin disk of cold-rolled steel, fixed to the 
housing (see Fig. 5). The magnetic fluxes of the magnet pass from the 
north pole of the magnet through the disk on both halves and back to 
the south pole, constituting a closed circuit. Except possibly for a small 
leakage, the unit does not act like a magnetic dipole with respect to the 
outside field. In the part of the disk near the poles of the magnet there 
is a relatively high and nonuniform magnetic field. Let the magnetic 
field in a magnetic domain i be H , , and let the induced magnetization 
in the same domain be /, , which is generally making an angle w with 
Hi . The magnitude of the retarding torque can be represented by 

f d = -SW,4hn. (39) 

i 

The minus sign means that the torque tends to oppose the relative 
displacement between the disk and the magnet. 

Provided that the spring constants lie in the ranges specified by (38), 
computer solutions have shown that there is a relatively small variation 
in damping time if the damping torques are in the ranges: 

T dl /I£? = 0.12 to 0.29 (roll) 

Tu/Itf = 0.16 to 0.45 (pitch) 

for the specified satellite. At 6000 nra, the numerical values of T d are 

T dl = 0.30 X 10" 4 to 0.72 X 10~ 4 ft-lb 

(41) 

T d2 = 0.40 X 10~ 4 to 1.12 X 10" 4 ft-lb. 

If the damping torques are much lower than the minimum values given 
above, the satellite will become earth-pointing only after a large number 
of orbits, as indicated by computer solutions. On the other hand, if the 
damping torques are much larger than the maximum values, the relative 
displacements become small and the satellite will keep tumbling for 
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many orbits. In the case of large angle motion the damping time is also 
found to be greater for values of t,i above the ranges given in (41). 

7.3 Hardware Development Program 

Damper units have been developed whose spring constants and 
damping torques fall in the range given by (38) and (41). These con- 
stants are suitable for the 6000-nm satellite previously described. How- 
ever, the mechanical design is such that the spring constants and 
damping torques may be adjusted for use with satellites at different 
altitudes (e.g., between (500 nm and 19,360 lira). 

7.3.1 Spring Design 

Successful torsion spring designs have evolved using both steel wires 
and flat beryllium-copper ribbons. The torsion springs must be under 
sufficient tension to prevent the lateral forces (due to gravity dif- 
ferentials, rotational motions, and environmental effects such as solar 
radiation, etc.) from deflecting the rotor laterally beyond the established 
clearance, thereby avoiding rubbing or sticking against the housing 
stops. The lateral forces have been calculated to be smaller than 10 -3 lb 
for the 6000-nm satellite previously described. An axial tension force of 
lb will be more than adequate to resist forces of this level. 

Two high-strength steel wires, each of 2-in. length and 0.008-in. diame- 
ter, will meet all of the specified requirements and provide a torsional 
spring constant of 0.36 X 10~ 3 ft-lb/rad. With a suitably designed end 
support for the springs, a number of torsional fatigue tests have shown 
that the springs are capable of withstanding in excess of \ million cycles 
at an amplitude of 60°. The static axial tension was 6 to 10 lb. This 
number of cycles is equivalent to 5 times the expected numbers of libra- 
tion periods in a 20-year useful life of the satellite. 

7.3.2 Damping Torque. Test Program 

The torque-displacement relationship (e.g., T*\ versus a for pitch 
displacement) has been measured for a damper unit at angular speeds 
between 0.5 X 10 -4 and 2 X 10~ 4 rad/sec (corresponding to a range of 
angular speeds of 0.18 9. to 0.73 9, at an altitude of 6000 nm). No de- 
pendence of damping torque on the angular speed has been observed 
in any of the tests, thereby verifying the assumption of velocity-inde- 
pendent magnetic hysteresis damping. Measurements have been made 
on a number of annular disks of various thicknesses made of cold-rolled 
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Fig. 13 — Rotational magnetic hysteresis loops. 

steel, annealed and miamiealed, and of V-Permendur. The maximum 
torque, T a , depends on the volume of the disk, the applied magnetic 
field and the degree of cold working of the material. A typical T d i-a 
curve measured on an unannealed cold-rolled steel (1010) is reproduced 
in Fig. 13. The energy dissipated per cycle is proportional to the area 
enclosed by the loops in the Tdi-a diagram. The slanted part of the 
curve extends over an angular displacement, 2a £b 8°, as shown. When 
the amplitude of the oscillation motion is less than a = 4°, minor loops 
as shown in Fig. 13 will be traced out. An appreciable loop area is still 
obtained even when a is as low as 1 °. 

Measurements have been made on the permeable disks after the disks 
have been irradiated by an electron flux of 10 17 /cm 2 , which is roughly 
equivalent to the highest electron radiation level anticipated within the 
Van Allen belt for a period of 30 years. The results indicate that electron 
radiation has very little effect on the damping torque. The effects of 
proton bombardment at a flux of 3 X 10 12 protons/cm 2 have also been 
found insignificant on the unannealed cold-rolled steel disks of 0.004- 
to 0.008-inch thickness without shielding. This flux is equivalent to 
the highest proton radiation level at (1000 nm for a period of 6 years. 

It has also been experimentally observed that the damping torque is 
relatively insensitive to wide temperature changes. The torque increases 
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only 10 per cent at — 40°F and decreases only 15 per cent at +250°F, 
from its value at room temperature. The temperature of the damper can 
be controlled within much closer limits in space by appropriate coatings, 
if desired. 

Vibration tests at a 20-g Level over a wide frequency band show that 
the damper will withstand launch conditions. 

VIII. ANALYSIS OF DISTURBANCES AND ERRORS 

In this section, we shall study the nature and magnitude of disturbing 
torques which produce forced librational motion of a gravitationally 
oriented satellite. In Sections 8.1 to 8.6, calculated values are given for 
each disturbing torque and its corresponding libration angle. 

In Section 8.7, the accumulative effects of all the disturbing torques 
are summarized. It will be seen that the satellite has been so designed 
that the gravitational torque dominates all disturbing torques at the 
altitudes of interest. 

It is obvious that the amplitude of steady-state librational motion 
should be kept to a minimum in order that maximum gain can be 
achieved from the earth-pointing antenna. For example, the theoretical 
(solid angle) gain at 0000 nm is 14.5 db with no allowance for librational 
motion. Allowance of a conservative tolerance of 10° on an antenna half 
angle, to accommodate 10° libration amplitude, results in a 3-db reduc- 
tion of theoretical antenna gain. However, it will be shown that the 
steady-state librational amplitude will be less than 10°. 

8.1 Solar Radiation Pressure 

An incident photon beam from the sun to a surface element will be 
partly absorbed, partly diffusely reflected and partly specularly reflected 
by the surface, resulting in an exertion of forces in directions normal and 
tangential to the surface element. These forces produce a net torque 
about the center of mass of the satellite. A detailed enumeration of the 
torques contributed by different surface elements on the two-body satel- 
lite shown in Fig. 3 indicates that a net maximum solar radiation torque 
of 0.5 X 10~ 4 ft-lb will act on the satellite. The magnitude of the gravita- 
tional torque at GOOO-nm altitude is 0.13 X 10~ 4 ft-lb per degree of angle, 
6, off the local vertical in the orbital plane for small libration angles 
(T 0mmx = |0 2 (/, - h) = 0.37 X 10" 3 ft-lb at 6 = 45°). Thus, statically 
the solar torque is balanced by the gravitational torque at 9 = 4°, 
when the sun is in its most unfavorable position. From computer solu- 
tions of the dynamics analysis in the pitch case, it is found that the 



2232 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1963 

satellite, which is provided with damping, will perform oscillations about 
the local vertical of a maximum amplitude not greater than 4°. Both 
the static and dynamic analyses neglected the rod deflections due to 
solar heating, and the accumulative effects due to solar torque and rod 
thermal bending will be discussed in the summary, Section 8.7. To be 
certain that there are no effects which would cause the libration ampli- 
tude to appreciably exceed 4°, it would be necessary to perform a three- 
dimensional analysis. 

The foregoing analysis was for a 6000-nm orbit employing the high- 
deck configuration. For significantly higher orbits the low-deck configura- 
tion would be preferred since it would experience less solar torque and a 
correspondingly smaller deviation from the local vertical. 

8.2 Residual Magnetic Dipok Moment 

The traveling-wave tube employed in a communications satellite 
such as the Telstar satellite contains two permanent magnets of equal 
size with the opposite poles placed against each other, thus constituting 
a quadrupole. Because of possible unequal strength of the two magnets 
and of inhomogeneous magnetic shielding outside of the traveling-wave 
tube, there would exist a net residual magnetic dipole moment in the 
satellite. Both the dipole and the quadrupole moments will interact 
with the geomagnetic field to produce torques. It can be shown that the 
torque produced by the quadrupole moment is only about 1 per cent of 
that produced by the residual dipole moment, when the moments of the 
two magnets are off by as little as 0.1 per cent. The magnetic moment 
of the Telstar satellite (produced mainly by the traveling- wave tube) 
was largely cancelled by the addition of compensating magnets. The 
residual dipole moment was 10 -6 weber-meter, the magnitude of which 
does not seem to have changed much after the satellite was launched into 
orbit. The use of two traveling-wave tubes, as might be needed in the 
commercial system, would not appreciably change the satellite's residual 
magnetic moment. 

Other magnetic dipole moments, which exist in the hysteresis damper 
units and the electric motors of the rod extension units, have been 
measured to be about 1.4 X 10 -6 weber-meter (a value obtained by 
adding all the moments scalarly). Therefore, a total magnetic dipole 
moment of 2.4 X 10~ 6 weber-meter may be expected in the satellite. 
This value may be reduced by "compensating" the motor dipoles and 
by further refinement of the cancellation techniques used on the Telstar 
satellite. Assuming a maximum geomagnetic field of 2.71 amp-turn/ 
meter (0.034 oersted) at 6000-nm altitude, we obtain a maximum torque 
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of about 0.65 X 10~ 5 newton-meter (0.48 X 10" 8 ft-lb). Upon balancing 
this torque against the gravitational torque in the manner indicated 
in the preceding Section 8.1, one finds a static libration angle of 0.4°. 

8.3 Orbital Eccentricity 

The eccentricity, e, of an elliptic orbit introduces a forcing torque on 
the satellite. If the eccentricity is not excessive, the satellite will settle 
down, as a result of damping, from an initial tumbling motion to a 
steady-state forced librational motion. In this case, the forcing torque 
due to eccentricity, 2e/£2 2 sin (fi/ + <p n ) (where / = /•> or = 7 B ), occurs 
only in the equations of pitch libration. In the case of viscous damping, 
the steady-state pitch librational angle has been found;"" however, 
because of the complexity of the resulting mathematical formulas, they 
will not be reproduced here. Since an analytical solution has not been 
obtained in the case of magnetic hysteresis damping, the steady-state 
libration angle of the earth-pointing body has not been evaluated 
exactly. However, it may be computed approximately by replacing the 
hysteresis damping by an equivalent viscous damping for the same 
energy dissipation per cycle (good only for small oscillations). In so 
doing, it is found that the libration amplitude 6 K, 3.6e radian, corre- 
sponding to the numerical data given in Section 2.2.4 and for spring 
constants and damping torques which fall in the range given in Section 
VII. A few computer solutions for the case of hysteresis damping indi- 
cate that 6 % 5e radian. Guided final-stage vehicles are believed to be 
capable of achieving orbit eccentricities below 0.005, which would result 
in libration amplitude of about 1.5°. 

s.4 Meteorite Impact 

Rased on Whipple's data," the meteorite flux rate to a spherical 
surface in the neighborhood of the earth can be shown to be 

$ = C/M meteorites per meter 2 per year (42) 

for meteorites of mass ^ M gram in the range of 10 -11 to 10" 1 gram. 
The constant C (in gram/meter 2 -year) is found to be C = 4.16 X 10~ 5 
according to Whipple and to be C = 20.8 X 10 _s according to Dubin. 22 
For meteorites hitting the deck body, which is at a distance L from the 
center of mass of the satellite, it can be shown that the expected number 
of meteorite collisions per year which result in satellite tumbling is 

■jr nvC 

N ' - 1 LA vm (43) 
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where v is taken to be the average speed of meteorites, A to be the 
average surface area of the deck body (neglecting the shadowing effect 
of the earth), /„ is the maximum moment of inertia of the composite 
satellite, p = (I x — I*)/Iy based on the composite satellite, and n is a 
factor determined by the momentum transfer (<1 for penetration, 
= 1 for completely inelastic impact, = 2 for perfectly elastic impact, 
and > 2 for hypervelocity impact when the material is blown backward 
out of a nearly hemispherical crater). Equation (43) is based on an 
analysis of the planar pitch motion of a single rigid body^atellite, which 
indicates that if an initial angular speed greater than V3pft is suddenly 
imparted to a satellite, which is already in line with the local vertical, 
the satellite will overcome a potential crest and turn over. Based on the 
result for pitch motion, the expected number of meteorite collision per 
year to give rise to angles of disturbance from the local vertical in the 
range from 0i to 2 ( ^90°) has been found to be 

% . nvC / 1 1 \ , . 

N > = Z LA V3^,n \sin ft ~ sin 6,) ' W 

Both (43) and (44) were derived in a simple manner by approximating 
the deck as a spherical body. A more exact result can be obtained if the 
meteorite flux rate is defined with respect to the projected area of a 
body. In this case the resulting expressions for i9\ and #2 are similar to 
(43) and (44), respectively, except that the coefficient (7r/4) LA is re- 
placed by complicated integrals involving the projected area element of 
various bodies and its distance from the mass center of the satellite. 

Numerical values of ft given in (43) and (44) calculated for the two- 
body satellite with C = 4.16 X 1(T B are tabulated in Table VIII, from 
which it is noted that the expected number of turnovers is 0.044 per year 
(or once in about 23 years). If C = 20.8 X 10" 6 is used, based on 
Dubin's 22 data, all values of ft in Table VIII should be multiplied by a 
factor of 5, and the expected turnover rate is 0.22 per year or once in 
about 4.5 years. These disturbances will be hysteretically damped 
down to a librational motion with amplitude of 5° in a reasonably short 
time. For example, computer solutions indicate that the pitch amplitude 
will be reduced from 45° to 5° in two to four orbital periods. In view of 
the uncertainty of the meteorite flux rate, all numerical values calcu- 
lated in this section are to be interpreted as order of magnitude estimates. 

8.5 Cocked Angle Due to Rod Deflections 

The extensible rods will be bent in a natural way and due to the ther- 
mal effects, as analyzed in Section V. Consequently, the axes of principal 
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moments of inertia of both bodies will deviate from their positions in 
the case of perfectly straight rods, and the center of mass of the com- 
posite satellite will not lie in the mast rod. As a result, the two bodies 
may not be perpendicular to each other (or the springs between them 
may not be in a neutral position) and the desired earth-pointing axis of 
the satellite may be off from the local vertical by a small cocked angle 
when the satellite is in a stable equilibrium position. Techniques have 
been developed for measuring rod straightness so that the natural rod 
bending in a gravity-free condition will not exceed a predetermined 
value. It has been found possible to select rods so that the rod bending 
will not exceed 1 foot for a rod length of (50 feet. This cocked angle 
could be evaluated if we could find the position vector of the hinge joint 
in the distorted configuration. A general error analysis has not been 
made, since it is not known a priori in what way the rods might be 
deflected. Based on the case of pitch libration, it is found that a deflec- 
tion of 1 ft of a (iO-ft long mast rod will cause a cocked angle of about 
0.7° in the case of the high-deck configuration. Assuming that the lateral 
tip deflections of both the mast and deck rods occur in the same direc- 
tion, the total cocked angle will be approximately 1.5°. For the case of 
the low-deck configuration, the cocked angle would be appreciably less. 
Rod bending can be caused by solar heating. As has been covered in 
Section 5.2, the deflection for a 50-ft long rod is expected to be about 2 
feet when the sun is perpendicular to the rod. There is a cumulative 
effect due to the mast and one pair of deck rods when the tips all bend 
away from the sun during certain periods of the year. The cumulative 
effects of these rods being bent produce a maximum cocked angle of 'A . 
This cocked angle can be reduced to less than 1° by silver-plating the 
rod exterior (the low absorptivity of silver, a ~ 0.1, would significantly 
reduce thermal rod bending due to a reduction of temperature dif- 
ferential across the rod cross section). As will be discussed in the sum- 
mary, Section 8.7, the effects of rod bending and solar torques are not 
additive. 

8.6 Miscellaneous Torques 

Torques due to self-gravity and eddy currents have been found to be 
much smaller than those discussed above. It can be shown that the self- 
gravity torque acting on one body due to the attraction of the other 
body is negligibly small compared to the gravitational torque at the 
altitude of interest. This is due to the fact that the sizes of the two bodies 
are not significantly different and that their masses are much smaller 
than that of the earth. Eddy-current losses induced in the conducting 
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materials of the satellite are small because of the low geomagnetic field 
at the altitudes of interest and because of the satellite's low angular 
speed. The torques due to self-gravity and eddy currents are less than 
10" 8 ft-lb. The torques produced by plasma effects due to the motion 
of the satellite in the Van Allen radiation belts are believed to be small. 
However, it is intended to make a detailed analysis of plasma effects. 

8.7 Summary 

In Table VIII, the maximum libration angles arc summarized for each 
individual disturbance for the high-deck configuration. By simply super- 
posing the individual effects of the various torques, the maximum 
steady-state libration amplitude is about 10° for the high-deck configura- 
tion. However, the various maximum individual effects cannot simply 
be added to obtain expected maximum libration amplitude. For example, 
the effect of solar torque and solar rod bending are not additive. The 
solar torque causes the deck assembly to rotate about the center of mass 
of the composite satellite in a direction away from the sun, whereas the 
rod bending due to solar heating causes an effect in the opposite direc- 
tion. A quantitative analysis is being made of the compensating effects 



Table VIII — Effects of Disturbances for a 6000-nm System 



Sources 



Maximum Magnitude 
(ft-lb) 



Approximate Librational 
Angle 



Pitch gravity torque: 1.3 X 10 _B per degree off the local vertical 



Solar radiation* 
Rod deflection* 

Natural bending 

Solar heating 
Orbital eccentricity 
Magnetic dipole moment 
Self-gravity and eddy current 



5 X icr s 

See text 



See text 

0.48 X 10-' 

Negligible 



1.5° 

3.0° 

1.5° for e - 0.005 

0.4° 

Negligible 



Effects 


* of Meteorite Impact 




Angle from the Local Vertical 


JV, Expected Number of 
Occurrences per Year 


Period of Occurrence in Years 


5°-15° 


0.336 


3 


15°-30° 


0.082 


24 


30°-50° 


0.031 


33 


50°-70° 


0.011 


94 


70°-90° 


0.003 


355 


5°-90° 


0.460 


2 


>90° (turnover or 


0.044 


23 


tumbling) 







* Computed for high-deck configuration. 
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of these two disturbances. Should it develop that the disturbances do 
not substantially compensate for each other, silver-plating the rods would 
reduce the effects of both disturbances to less than 5°, rather than 7°, 
which is the sum of the two disturbance angles. The effects of the other 
disturbances — natural rod bending, orbit eccentricity, magnetic dipole 
moment — would not be added to give a total angle of 3.4°. Hence for 
the high-deck configuration, the final librational angle is expected to be 
well under 10°. For the low-deck configuration, the librational angle 
would be expected to be somewhat smaller than that for the high-deck 
configuration. 

IX. CONCLUSIONS 

The theoretical feasibility of the proposed PGAC system has been 
amply demonstrated by more than one hundred computer runs based 
on the dynamics analysis 4 of the ideal two-rigid-body system. Computer 
simulations made with a wide variety of initial conditions showed that 
the system stopped tumbling and then, within about 7 orbital periods, 
settled down to a state of small oscillations about an earth-pointing 
direction. It has also been indicated that disturbing influences, such as 
solar radiation pressure and orbital eccentricity, produce oscillations 
of less than 10° for a (iOOO-nm orbit. 

The rods have been shown to possess adequate rigidity, to be fully 
capable of withstanding the loads imposed during the extension phase, 
and not to undergo excessive bending due to solar heating. 

On the basis of comprehensive studies and tests, it is believed that 
the PGAC system described in this paper is fully capable of meeting 
all its design objectives, including compatibility with multiple launch 
procedures, and that it will provide a significant advance in communica- 
tions satellites practice. 
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